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We study the kinetics of the quasiparticle capture and emission process in a small superconduct- 
ing island (Cooper-pair box) connected by a tunnel junction to a massive superconducting lead. At 
low temperatures, the charge on the box fluctuates between two states, even and odd in the number 
of electrons. Assuming that the odd-electron state has the lowest energy, we evaluate the distribu- 
tion of lifetimes of the even- and odd-electron states of the Cooper-pair box. The lifetime in the 
even-electron state is an exponentially distributed random variable corresponding to a homogenous 
Poisson process of "poisoning" the island with a quasiparticle. The distribution of lifetimes of the 
odd-electron state may deviate from the exponential one. The deviations come from two sources - 
the peculiarity of the quasiparticle density of states in a superconductor and the possibility of quasi- 
particle energy relaxation via phonon emission. In addition to the lifetime distribution, we also find 
spectral density of charge fluctuations generated by capture and emission processes. The complex 
statistics of the quasiparticle dwell times in the Cooper-pair box may result in strong deviations of 
the noise spectrum from the Lorentzian form. 



I. INTRODUCTION 

Properties of a mesoscopic superconducting circuit 
may depend crucially on the presence of quasiparticles in 
its elements. The operation of a superconducting charge 
qubit, for example, requires two-electron periodicity of 
its charge staie^i^^^^^^. This periodicity may be in- 
terrupted by the entrance of an unpaired electron into the 
Cooper-pair box (CPB) serving as an active element of a 
qubit. The quasiparticle changes the charge state of CPB 
from even to odd, and lowers the charging energy. This 
trapping phenomenon, commonly referred to as "quasi- 
particle poisoning" , is well-known from the studies of the 
charge parity effect in superconductors^^^^ . Quasiparticle 
poisoning contributes to the phase relaxation in super- 
conducting qubits^^. For a typical CPB size and tunnel 
conductances of the order of unit quantum, the quasipar- 
ticle dwelling times are of the order of a few /is. This time 
scale is at the edge of accessibility for the modern exper- 
iments^. Individual quasiparticle tunneling events were 
resolved and the statistics of quasiparticle entrances and 
exits from the CPB box were investigated in Refs. 00] • 
The observed statistics of entrances was well described by 
a standard Poissonian process^^^. For the quasiparticle 
exits, the results are less clear. In many cases, it may be 
well described by the Poissonian statistics^i^. However, 
there are indications of deviation from that simple law 
for some samples 

In this paper, we develop a kinetic theory of quasipar- 
ticle poisoning. We find the distribution of times A/odd(^) 
and A^even(^) the CPB dwells, respectively, in odd- and 
even-electron states. We also find the spectrum of charge 
noise produced by the poisoning processes. 

The conventional Poissonian statistics of the quasipar- 
ticle exits would yield an exponential distribution for 
odd-electron lifetime in the box. We see two reasons for 
the distribution function A^odd(^) to deviate from that 
simple form. The first one is related to the thermaliza- 
tion of a quasiparticle within the CPB. If the rates of 



energy relaxation and of tunneling out for a quasipar- 
ticle in CPB are of the same order, then two different 
time scales control the short-time and long-time parts 
of the distribution function A/odd (^)- The shorter time 
scale is defined by the escape rate Tout of the unequili- 
brated quasiparticle from the CPB. The longer time scale 
is defined by the rate of activation of equilibrated quasi- 
particle to an energy level allowing an escape from CPB. 
The second reason for the deviations from the exponen- 
tial distribution controlled by a single rate, comes from 
the singular energy dependence of the quasiparticle den- 
sity of states in a superconductor. Because of it, the 
tunneling-out rate depends strongly on the quasiparticle 
energy. Thus, even in the absence of thermalization the 
quasiparticle escapes from CPB cannot be described by 
an exponential distribution. 

The conventional Poissonian statistics for both en- 
trances and exits of the quasiparticle would lead to a 
Lorentzian spectral density Sq{u;) of CPB charge fluctu- 
ations^^. The interplay of tunneling and relaxation rates 
may result in deviations from the Lorentzian function. In 
the case of slow quasiparticle thermalization rate com- 
pared to the quasiparticle tunneling-out rate Tout, the 
function Sq{uj) roughly can be viewed as a superposition 
of two Lorentzians. The width of the narrower one is 
controlled by the processes involving quasiparticle ther- 
malization and activation by phonons, while the width of 
the broader one is of the order of the escape rate Tout- 

The paper is organized as follows. We begin in Sec. [Ill 
with the qualitative derivation and discussion of the main 
results. In the next sections piIIIIVp we derive and solve 
the microscopic master equations for the kinetics of the 
quasiparticle capture and emission, and calculate the life- 
time distribution functions in the even- and odd-charge 
states. In Sec.|V]we calculate charge noise spectral den- 
sity Sq{uj) for the Cooper-pair box. In SecEBwe sum- 
marize the main results. Some technical details are rele- 
gated to the Appendix. 



FIG. 1: (color online). Schematic picture of the Cooper-pair 
box qubit. The left-hand side superconducting mesoscopic is- 
land is the Cooper-pair box connected via a tunable Josephson 
junction to the large superconducting lead (right-hand side). 
Gate bias is applied through the capacitance Cg. The junc- 
tion is characterized by the dimensionless conductance . 



II. QUALITATIVE CONSIDERATIONS AND 
MAIN RESULTS. 

A. Relevant time scales. 

Dynamics of the Cooper-pair box coupled to the su- 
perconducting lead through the Josephson junction, see 
Fig. [H is described by the Hamiltonian 

H = ^ H^cs + ^Bcs + • (1) 

Here i^BCS ^bcs BCS Hamiltonians for box and 
the lead; = EdQ/e-Ng)^, with Ec, Ng, and Q be- 
ing the charging energy, dimensionless gate voltage, and 
charge of the CPB, respectively. The tunneling Hamilto- 
nian Hrj. is defined in the conventional way 

Ht=Y. {^"A.cp,, + H-c.) , (2) 

where tkp is the tunneling matrix element, Ck^o- and Cp^cr 
are the annihilation operators for an electron in the state 
|/c,cr) in the CPB and state |p, cr) in the superconduct- 
ing lead, respectively. Here superconducting gap energy 
is the largest energy scale, A > Ec > Ej T. In or- 
der to distinguish between Cooper pair and quasiparticle 
tunneling, we present the Hamiltonian ([T]) in the form^^ 

H = Ho^V, SindV = H^-H^. (3) 

Here Hq = ^ H^^^^ + H^^^^ + Hj, and Hj is the 
Hamiltonian describing Josephson tunneling 

H, = \N) (iV|JT^-_l_ifJiV+2) (iV+2|+H.c. 

The matrix element (A^| e-Hq \^ + 2) propor- 
tional to the Josephson energy Ej. (Here Ej is given 
by the Ambegaokar-Baratoff relation.) The perturbation 
Hamiltonian V defined in Eq. ([3]) is suitable for calcula- 
tion of the quasiparticle tunneling rate. 

Energy of the system as a function of the gate voltage 
is shown in Fig. [2l Ai Ng = 1 the electrostatic energy 
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FIG. 2: (color online). Energy of the Cooper-pair box as a 
function of dimensionless gate voltage Ng in units of e. Solid 
(red) line corresponds to even-charge state of the box, dashed 
(blue) line corresponds to the odd-charge state of the box. 
The trap depth 8E is the ground state energy difference be- 
tween the even-charge state (no quasiparticles in the CPB), 
and odd-charge state (an unpaired electron in the CPB) at 
Ng — \. (We assume here equal gap energies in the box and 
the lead, = = A.) 

of the system is minimized when unpaired electron re- 
sides in the CPB. Thus, at Ng = 1 the CPB is a trap for 
a quasiparticle. The trap depth ^E is equal the ground 
state energy difference between the even-charge state (no 
quasiparticles in the CPB) and odd-charge state (an un- 
paired electron in the CPB). For equal gap energies in the 
box and the lead, = A5 = A, the trap is formed due 
to the Coulomb blockade effect. In the case Ec ^ 
one has 

5E^E,-^:^ T, (4) 

and only two lowest charge states are important, see 
Fig. O Also, we assume that there is at most one quasi- 
particle in the box in the odd state^^. 

The transition probability per unit time between odd 
and even-charge states W{p^ k) can be obtained using the 
Fermi golden rule {h= 1), 

W{p, k) = 2^\ (p, e| V |o, k) fS{Ep^SE-Ek). (5) 

Here the state |e,p) corresponds to even-charge state of 
the box and the quasiparticle in the state \p ) in the reser- 
voir; the state |o, k) corresponds to the odd-charge state 
of the box and quasiparticle in the state \k) within the 
box. The quasiparticle energies in the CPB and lead 

Ef./p are defined as Ej./p = ^^^^/^ + • Matrix ele- 
ments (p, e| V|o, /c) can be calculated using the Bogoli- 
ubov transformationi^iii. Taking into account the rela- 
tion between tunneling matrix elements and the normal- 
state junction conductance the expression for W{p^ k) can 
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FIG. 4: The dependence of the escape rate Tout (Ek) on energy 



FIG. 3: (color onhne). Schematic picture of the CPB-lead sys- 
tem showing ahowed transitions for the quasiparticle injected 
into the excited state of the box. At Ng = l the Cooper-pair 
box is a trap for quasiparticle. 



be written as 



Wip, k) = 



9t 



Stt 



l + ^-^^-^5{E, + 5E-Eu) (6) 



with 5i)ii being mean level spacing in the box/lead, and 
Qrj. being dimensionless conductance of the junction. 

Using the transition rate (|6]), one can calculate the 
level width of the state |o, k) with respect to quasiparticle 
tunneling through the junction to the lead, 



Tont{Ek) = Y.W{p,k) 

p 

Air {Ek-SE)Ek 

The Heaviside function Q{x) appears in Eq. ([7j) because 
there are no states to tunnel into for a quasiparticle with 
energy lower than the threshold energy £^thd, see Fig. [3l 



(7) 



v{Ek-dE)e{Ek-Eti,i). 



^thd 



A + 5E. 



(8) 



The quasiparticle density of states v{Ek) (in units of the 
normal density of states at the Fermi level) is given by 



(9) 



Due to the square-root singularity here, the rate rout(^/c) 
has square-root divergence at Ek = £^thd, see Fig. (j4]). 

The quasiparticle may enter and subsequently leave the 
Cooper-pair box without changing its energy. For such 
elastic process, the excess energy of the exiting quasi- 
particle is equal to its initial energy, and is of the order 
of the temperature, i.e., Ek — Et\^^r^T. Therefore, the 
corresponding escape rate is 



Here for brevity we denote fiT) = u{Ek = A-\-T). For the 
system with < 1, volume of the CPB Vb < tem- 
perature T ~ 50mK and dE ~ 0.5K, the typical escape 
time r~J^ is of the order of a /is. 

To find the average rate Fin of quasiparticle tunneling 
from the lead to the CPB, we integrate the transition 
probability per unit time (|6|) with the distribution func- 
tion f{Ep) of quasiparticles in the lead. 



J2w{k,p)f{Ep 



(11) 



p,k 



Upon elastically tunneling into the excited state in the 
CPB the quasiparticle can relax to the bottom of the 
trap, see Fig. [3l For that, the quasiparticle needs to 
give away energy ^ SE. At low temperatures the domi- 
nant mechanism of quasiparticle energy relaxation is due 
to electron-phonon inelastic scattering rate l/r{Ek). At 
low temperature quasiparticles are tunneling into the box 
through the energy levels just above the threshold energy 
Ek ~ ^thd, see Eq. (|8]). Assuming 5E <C A, the typical 
quasiparticle relaxation time r is given by^^ 



r = r(E/e~£'thd)^ To 



SE_V 



(12) 



Here tq is a characteristic parameter defining the average 
electron-phonon scattering rate at T = Tc with Tc being 
superconducting transition temperature. In aluminum, a 
typical material used for CPB, tq ^ 0.1 - 0.5 
As one can see from Eq. (fT2]) , the quasiparticle relaxation 
rate is a strong function of the trap depth SE. There- 
fore, depending on 6E there are two kinds of traps - "shal- 
low" traps corresponding to rFout ^ 1^ and "deep" traps 
with rFout <^ 1- (Note, for shallow traps we still assume 
SE ^ T.) The important quantity characterizing the 
traps is the probability Ptr for a quasiparticle to relax to 
the bottom of the trap before an escape. 



(10) 



(13) 
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B. Lifetime distribution function. 

Experimentally observable quantity^^^, which reveals 
the kinetics of quasiparticle trapping, is the lifetime dis- 
tribution function A/odd (^) of odd-charge states of the 
CPB. The distribution of lifetimes A^odd(^) depends on 
the internal dynamics of the quasiparticle in the CPB, 
i.e the ratio of rPout- 

We start with the discussion of the long-time asymp- 
tote of the lifetime distribution function. At t > r the 
dwell-time distribution A^odd(^) is governed by phonon- 
assisted activation of the thermalized quasiparticle in the 
trap. The phonon adsorption processes are statistically 
independent from each other. Hence, the lifetime distri- 
bution exponentially decays with time 

Nodd{t) cx exp(-7t) (14) 

with the rate 

This expression can be understood as follows. The rate of 
thermal activation of the quasiparticle from the bottom 
of the trap to the threshold energy is 77^^^ exp (— ■^) , 
for brevity we define v{5E) = p{Ek = ^'thd)- The addi- 
tional factor Tj{5E) /v{T) here comes from the difference 
of the quasiparticle density of states at the bottom of the 
trap i^{T) and at the threshold energy v{5E). The last 
term (1 — Ptr) in Eq. (p!5|) corresponds to the probabil- 
ity of the quasiparticle escape to the lead upon activa- 
tion. Equation ([15]) allows us to consider limiting cases 
of rEout < 1 and rEout > E 

In the case of "deep" traps (rEout ^ 1) most quasipar- 
ticles upon entering the excited state in the box quickly 
thermalize . Therefore, the main contribution to lifetime 
distribution function comes from phonon-assisted escapes 
described by Eq. ([Hj), see Eig. ([5]). The activation escape 
rate of Eq. ([15]) in this limit equals 

^ u(5E) ( SE\ , , 

since 1 - Ptr ^ Eout^, see Eq. (fT3|) . 

In the opposite limit rEout ^ 1, ^-e. "shallow" traps, 
the probability for a quasiparticle to relax to the bottom 
of the trap is small Ptr 1- Therefore, upon elasti- 
cally tunneling into the excited state in the CPB the 
quasiparticles will predominantly return to the reservoir 
unequilibrated. Nevertheless, there is a small fraction of 
quasiparticles (~ 1/rEout) that do relax to the bottom of 
the trap, and stay in the box much longer than unequi- 
librated ones. Thus, at t > r the dwell-time distribution 
function 7Vodd(^) has an exponentially decaying tail ([H]), 
see Eig. ([5]), with phonon-activated escape rate 

lu{SE) f SE\ 
7s--^exp^--j. (17) 



At t ~ r the typical value of the lifetime distribution 
function is Nodd{t^r) ^ 7s /rEout- 

At short times, t <^ the lifetime distribution 
function Nodd{t) describes the kinetics of unequili- 
brated quasiparticles. Quasiparticles tunnel into the box 
through the energy levels E^ = £^thd + ^ (here £ > 0), 
and predominantly reside there until the escape with the 
rates Eout(^)- For a given energy level s the lifetime dis- 
tribution is exponential 

Nodd{£.t) (X exp(-Eout(s:)^). (18) 

Note that upon entering into the CPB from the reservoir 
the quasiparticles can populate different levels within the 
energy strip ~ T, see Eq. (|36|) . Therefore, experimentally 
observable quantity A^odd(^), obtained by the statistical 
averaging over a large number of the tunneling events, is 
given by 

Noddit) ^ deexp (^-^ - Eout(e)^) • (19) 

Taking into account the singularity of Eout(^) at small 
energies Eout(^) ^ s~^l^ ^ we find that A/odd (^) deviates 
from the simple exponential distribution [see Eig. ([5])], 

iVodd(t) oc exp 1^-3 {^^^ j (20) 
at times t > 1/Eout- See also Sec. ([IV]) for more details. 

C. Charge Noise Power Spectrum. 

Anomalies in the lifetime distribution, see Eig. ([5]), 
should also lead to a specific spectrum of charge fluc- 
tuations. We define the spectral density of charge fluc- 
tuations Sq{u;) in the Cooper-pair box as 

/CO 
dte"^' [{SQit)SQiO)) + {SQ{0)SQm 
-OC 

(21) 

with 5Q{t) = Q{t) — (Q). The variance of the fluctuations 
of charge Q in the CPB, 

= I -^Sq{lo), (22) 

is a thermodynamic, not a kinetic, quantity, and is known 
from statistical mechanics. The kinetics of the system is 
reflected in the dependence of the noise spectrum ([2T]) on 
the frequency co. 

In the limit of fast relaxation rEout 1 the es- 
capes from the CPB are given by one time scale ([T6]) . 
The quasiparticle entrances into and exits from the CPB 
are random, and can be described by Poisson processes. 
Thus, Sq{u;) is given by the Lorentzian function corre- 
sponding to random telegraph noise—, 

SQ{ij) ^ 4e^a-odd(l - ^odd)7 I 1 • (2^) 
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FIG. 5: (color online), a) Schematic picture of the life- 
time distribution function for "deep" traps (rFout ^1). b) 
Schematic picture of the lifetime distribution function for 
"shallow" traps (rFout ^ !)• Inset: Deviations of Noddif) 
from exponential distribution at short times. 



The second (quasi) Lorentzian function Sq{(jo) is associ- 
ated with fast charge fluctuations reflecting the kinetics 
of unequilibrated quasiparticles. Assuming uo ^ Fout ^ 

(2) 

Fin the asymptote of Sq{uj) is 

S^^Hco) ^ e^^exp{-SE/T) (^)\ (26) 
^ Tout V ^ / 

(2) 

The width of Sq (co) is determined by the typical es- 
cape rate of unequilibrated quasiparticles from the box 
Fout defined in Eq. (p!Q|) . Similar to the lifetime distribu- 
tion, see Fig.[5l we predict deviations of 5'g^(a;) from the 
Lorentzian function at ~ Fout due to the peculiarity of 
the quasiparticle density of states. 

The high-frequency tail of Sq{u;) is provided by 

Eq. ([26|) . However, the contribution of Sq {co) to the 

sum rule ([22|) is much smaller than that from Sq\uj). 
In other words, the main contribution to the noise power 
comes from slow fluctuations. It resembles the case of 
the current noise in superconducting detectors^^. 

In the rest of the paper, we provide detailed derivation 
of the results discussed qualitatively in this section. 

III. LIFETIME DISTRIBUTION OF THE 
EVEN-CHARGE STATE. 



Here ^odd is an equilibrium average occupation of the 
odd-charge state in the CPB (0 < ^odd ^ 1), see 
Eq. ^Ej for details. At low temperature (T <C SE) 
the box is predominantly in the odd-charge state, i.e 
(1 — ^odd) ex.p{—SE/T). The rate of activated quasi- 
particle escape processes has the same small exponent, 
therefore the width of the Lorentzian ([23|) is mainly given 
by the transitions from even to odd-electron state, 

-!-=»r,.,. 

Teff 



(24) 



See Eq. ([83]) for the full result. 

In the limit of slow relaxation (rFout ^ 1) the charge 
noise power spectrum Sq{u;) deviates significantly from 
the Lorentzian. These deviations stem from the fact that 
a quasiparticle may escape from the box before or after 
the equilibration, which results in two characteristic time 
scales for the escapes^^, see Fig. [5l Consequently, the 
function Sq{u;) can be roughly viewed as a superposition 
of two Lorentzians, and is similar to carrier concentra- 
tion fluctuations in semiconductors due to trapping^^. 
The "narrow" Lorentzian describes the dynamics of slow 
fluctuations due to phonon-assisted trapping of quasipar- 
ticles 



2^odd(l - CToddj'Teff 1 1 Lij 



Teff r Fin + Fo 



-.(25) 



(Cjreff)2 + 1 

The width here is determined by the probability of 
quasiparticle trapping per unit time. (Like above, we as- 
sume here T <^ 6E and neglect activated escape rate.) 



Let us assume that the system switched to the even 
state at t = 0, and introduce the probability density 
A^even(^/c7 ^) for ^ quasiparticlc to enter the CPB for the 
first time through the state E}^. Then, the probability 
density for the CPB to reside in the even state until time 
t is 



Neven{t) = ^ A^even , ^) • 



(27) 



A'evenl^fc, ^) is givcu by the conditional probability of 
quasiparticle entering the CPB into an empty state Ek 
during the interval (t, t + dt) times the probability that 
any quasiparticle has not entered into any state in the 
CPB during the preceding interval (0,t), 

N,^,n{Ek,t)dt = Y,W{k,p)f{Ep) 

p 

X (^^-Y^l dt'N.^UEk' , t')^ 428) 

Summing Eq. ([28]) over states k and solving for A'even(^) 
one finds 



Neven{t) = Fin exp (-Pint) , 



(29) 



which corresponds to a homogenous Poisson process. The 
quasiparticle tunneling rate from the lead to the CPB Fin 
is given by Eq. (pT]) . 

Recent experiments by Aumentado et. al. [HQ] indi- 
cate that the density of quasiparticles in the lead 



6 



exceeds the equilibrium one at the temperature of the 
cryostat. The origin of nonequihbrium quasiparticles is 
not clear, but it is plausible to assume that quasiparticle 
distribution function in the lead f{Ep) is given by the 
Boltzman function 



f{Ep)= exp 



En 



T 



(30) 



with some effective chemical potential and temperature, 
fii and T, respectively. The chemical potential /j^i is re- 
lated to the quasiparticle density by the equation 



(31) 



Here Vi is the volume of the lead. We consider the density 
of quasiparticles n^p and their effective temperature as 
input parameters here, which can be estimated from the 
experimental data^^^^^. Taking into account Eq. ([3Q|) we 
can evaluate the right-hand side of Eq. (pT]) to obtain 



qp 



v{8E) 



SE 



A^SE' 



(32) 



Here up is the normal density of states at the Fermi level. 
The average waiting time in the even-charge state is 



POO 

{reven)= / N^.^n{t)tdt = Tr\ 

Jo 



(33) 



This result is expected for the conventional Poisson pro- 
cess. 



IV. LIFETIME DISTRIBUTION OF THE 
ODD-CHARGE STATE. 

A. Master equation for survival probability. 

The distribution of dwell times for odd-charge state is 
more complicated than for even-charge state due to the 
internal dynamics of the quasiparticle in the CPB. Upon 
tunneling elastically into the box the quasiparticle enters 
into the excited state with typical excess energy SE above 
the gap in the island. The dwell time of the quasiparticle 
in the box depends whether upon tunneling into the ex- 
cited state it relaxes to the bottom of the trap or tunnels 
out un-equilibrated, see Fig. 2. In order to describe the 
physics of quasiparticle tunneling we develop a formal- 
ism similar to the rate equations theory. We will include 
electron-phonon collision integral into our equations to 
account for the internal dynamics of the quasiparticle in- 
side the CPB. The experimentally accessible quantity is 
the probability density A^odd(^) of leaving an odd state 
in the time interval (t, t -h dt) assuming that the quasi- 
particle resided continuously in the box during the time 
interval (0, t). The object convenient for evaluation is the 
conditional probability 6'odd(^) (or survival probability) 



for a quasiparticle to occupy a given level, under the con- 
dition that the unpaired electron continuously resided in 
the box over the time interval (0,t). The lifetime distri- 
bution A^odd(^) can be easily obtained from 5'odd(^), 

^odd(i) = |(1 - SoMt)) = (34) 

Probability 5'odd(^) is simply related to the conditional 
probability S{Ek, t) for a quasiparticle to occupy level E^ 
at the moment t in the box assuming that a quasiparticle 
entered CPB at t = and resided continuously in the box 
during the time interval (0,t), 



5'odd(^) = So(i(i{Ek,t). 



(35) 



We assume that in the initial moment of time the 
quasiparticle has just entered the state Ek in the box. 
Therefore, the initial probability S'o(£^/e,0) of the occu- 
pation of the level E^ in the box is determined by the 
tunneling rate into the state 



5o(i?fe,o) = -^^vr(p,fc)/(i;p). 



(36) 



The normalization of So{Ek^O) is chosen to satisfy 
5'odd(0) = Efe^o(£^/c,0) = 1. According to Eq. (|36]) the 
initial conditional probability So{Ek,0) is zero below the 
threshold energy Ek < ^thd, and is proportional to Gibbs 
factor above the threshold Ek > ^'thd- This reflects out- 
of-equilibrium quasiparticle distribution at t = 0. 

The conditional probability Sodd{Ek, t) consistent with 
initial conditions (|36|) satisfies the following master equa- 
tion: 



So{Ek,t)-^Tout{Ek)So{Ek,t) = 



Sodd{Ek,t)- 



(37) 



The second term on the left-hand side corresponds to 
the loss from the state Ek due to the tunneling through 
the junction to the lead with the rate rout(^/c) of Eq. ([7j). 
Note that unlike in the theory of the rate equations there 
is no "gain" term in Eq. (|37j) . This is due to the condi- 
tion that the box is occupied at t = and remains occu- 
pied continuously until some time t. The right-hand side 
of Eq. (|37j) corresponds to the electron-phonon collision 
integral in the relaxation time approximation with r of 
Eq. ([12]) and 

s:uEk.t)=Sodd{t)-p'odd{Ek)- 

Note that Eq. ([37|) is nonlocal in Ek due to the depen- 
dence of the collision integral on 6'odd(^) (see Eq. (j35j)). 
The collision integral in Eq. (|37j) describes the phonon- 
induced relaxation of the trapped quasiparticle to an 
equilibrium. 



exp j-Ek/T) 

Zodd 



(38) 
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Here T is the quasiparticle temperature in the box. 
(For simphcity, we assume that the effective quasipar- 
ticle temperature in the lead is the same as in the box, 
Ti = Tb =T.) The normalization factor Zodd at T <C T* 
is given by 



Zodd 




exp 



(39) 



B. General solution for Sodd{t). 

Using Laplace transform, 

So{Ek,s) = / dtSo{Ek,t)e-'\ (40) 

we reduce the differential equation (|37j) supplied with the 
initial conditions (|36|) to an algebraic one 



(41) 

Equation (|4T]) can be solved for 6'o(£^/c,s). Then, by 
summing that solution over momenta k and utilizing 
Eqs. (|35|) and ([36|) we find the survival probability 

Sodd{s) 



Sodd{s) 



B{s) 



1-A{s)' 

Here functions B{s) and A{s) are defined as 



(42) 



Bis) 



Tin ^ 5+l/r + rout(^0 



(43) 



^ ^ r ^ s + 1/r 



/r + rout(^fc)' 



At T ^5 one can take the thermodynamic limit and 
replace the sums with the integrals in Eq. (|43|) . Further 
simplification of the denominator in Eq. (|42|) is possible if 
one splits the integral in A{s) into the intervals (A, £^thd), 
where Tout{Ek) = 0, and (£;thd,oo). Then, Eq. (|42]) be- 
comes 



Sodd{s) 



1 

5+ - 



B{s) 



r J s^X{s) 
with the functions B{s) and X{s) defined as 



(44) 



B{s) 
X{s) 



L 



in JBthd 



dEk f{Ek-5E)rodd{Ek) 



2 r°° 



dEk . ^ p^^^{Ek)Todd{Ek) 



s^l/T^Todd{Ek)' 



(45) 




Re(s) 



FIG. 6: (color online). Contour of integration (red line) cho- 
sen to calculate inverse Laplace transform Eq. (|46p . Points of 
non-analytic behavior of cr^^ (s) are shown. Poles at si, S2, 
and a cut s G (— oo, — Smin). 



The inverse Laplace transform is given by 

ne-\-ioo 



So{t)-- 



1 



dsSodd{s)e^ 



(46) 



where e is chosen in such way that Sodd{s) is analytic 
at Re [5] > e. The integral (|46|) can be calculated using 
complex variable calculus by closing the contour of in- 
tegration as shown in Fig. [6] and analyzing the enclosed 
points of non-analytic behavior of Sodd{s)- In general, 
the singularities of 6'odd('5) consist of two poles and a 
cut. The latter is due to the singularities of the function 
B{s) causing 6'odd(<5) to be non- analytic along the cut 
s G (—00, — 5min), where 



Smin = - + min [ Fodd {Ek)]- 

T 



(47) 



The plot of Fodd(^/c) is shown in Fig. (|4|). The function 
^oddlEk) has a minimum at = £;thd + 5E/2. (For 

the estimate of the minimum we assumed 5E <C A.) In 
addition to the cut, 6'odd(<5) has 2 poles. The poles si 
and 52 are the solutions of the following equation in the 
region of analyticity of B{s)\ 



s^X{s) = 0. 



(48) 



We now analyze the singularities 5'odd(<§) and find their 
contribution to the integral (|46|) . 

The contribution from the cut to Eq. (|46|) corresponds 
to the kinetics of unequilibrated quasiparticles. Formally 
it comes from the non-analyticity of 5'odd(<§) due to the 
singularities of the function B{s) itself. The proper con- 
tribution to Eq. (|46|) can be calculated by integrating 
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along the contour enclosing the cut, 
-1 

hut = ^ / ^^^^^ {Sodd{s^ie)-Sodd{s-ie)) . 

^ Smin 

(49) 

At low temperature T <C SE^ the discontinuity of the 
imaginary part of 6'odd(<5) at the cut is 

SoM{s + ie)-Soid{s-ie) = 2i (^s + ^^ I]^^^^±!£l. (50) 

Substitution of this expression into Eq. (|49|) yields 

/cut = ^ [ ^ iy{Ek)f{Ek-SE)Todd{Ek) 

^ TT-F^^^^P ^odd{Ek)t]. (51) 

To simplify the above expression we introduce the dimen- 
sionless variable 

Ek — Ethd 



(52) 



and write the integral in /cut in terms of 



1 

^/7^Toutl'{oE) Jq 1 
X exp {-z - rodd(^)^ - t/r) . 



rT oddjz) 
^ rT odd{z) 

(53) 

Here and thereafter rodd(^) and ^{z) are given by 
Eqs. ([7j) and ([9]), respectively, with Ek = £^thd + Tz. 

We now analyze the contribution to Eq. (|46|) from the 
poles. The pole at si may be found from the iterative 
solution of Eq. (|48|) at small s (s <C Smin), 

si « -X(s = 0), (54) 

with X{s) given by Eq. (|45|) . The contribution from the 
pole at si can be easily calculated using residue calculus 
yielding 

h =r(0)exp[-X(0)t]. (55) 

Equation ([55|) describes the kinetics of thermalized quasi- 
particles. At low temperature X(0) (x exp(— (5£^/T), 
which justifies the approximation used in Eq. ([54|) . see 
also next section. The function Y{0) depends on rPout, 
and is approximately given by 



r(o) 



1 



dz 



exp {—z) 



Jo V^ + 'Tfout 



(56) 



Here we used small-z asymptote {z <^ ^) for the escape 
rate, 



rout(z) 



(57) 



The second pole S2 is given by the solution of Eq. ([48|) 
at large s. At small temperature T <^ dE one can show 
that the contribution of the second pole 52 to Eq. (|^6|) 
is small, and thus can be neglected. (For details, see the 
Appendix in Ref. [11] ) 



C. Results and Discussions. 

Combining all relevant contributions to the inverse 
Laplace transform, Eqs. ([53|) and ([55|) . we obtain the 
solution for the survival probability 



Soddit) = y(0)exp(-7t)+F(i). 



(58) 



The first term here corresponds to the kinetics of the 
quasiparticle that relaxed to the bottom of the trap. The 
thermally activated decay rate 7, found with the help of 
Eqs. ([511) and ([52]), is 



7 = - 



1 u{SE) 



exp 



SE 



dz 



e-^/r 



i/r + rout/v^y 



(59) 



The integral in Eq. ([59]) refiects the probability for a 
quasiparticle to relax to the bottom of the trap [cf. 
Eq. (fTS])]. The second term in Eq. ([58]) describes the 
kinetics of unequilibrated quasiparticles with F{t) given 

by 

Trodd(^) 

+ TTodd{z) 

(60) 



X exp -z - trodd(^) 



In the next paragraphs we will analyze 6'odd(^) for fast 
and slow relaxation limits. 

In the fast relaxation limit rPout 1 ("deep" 

trap), the leading contribution to the survival probability 
5'odd(^) comes from the first term in Eq. ([58]) . the second 
term in Eq. ([58]) is proportional to rPout , and can be ne- 
glected. Consequently, the survival probability is given 
by 

Sodd{t) - exp (-7ft). (61) 
Using Eq. we find the lifetime distribution function 

iVodd(t) =7fexp(-7ft), (62) 



cf. Eqs. (|T4|) and (|T6l) . As discussed qualitatively 
in Sec. [Til in the fast relaxation limit the majority of 
quasiparticles entering the CPB into the excited state 
Ek ~ ^thd relax to the bottom of the trap and stay in 
the box until they are thermally excited out of the trap 
by phonons with the rate 7f of Eq. ([TB]) . Finally, using 
Eq. ([61]) we find the average lifetime of the odd-charge 
state 



(^odd) 



JO 



Sodd{t)dt = l/7f 



(63) 



In the opposite limit of the "shallow" trap, rPout ^ 1, 
the majority of quasiparticles tunnel out unequilibrated 
to the lead (Ptr ~ l/rFout)- The expression for the sur- 
vival probability ([58l) in this limit becomes 



Sodd{t)=F{t) 



1 



V^rEo 



■exp(-7st). 



(64) 
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Note that in addition to the first term describing the ki- 
netics of unequihbrated quasiparticles the survival prob- 
abihty has a tail corresponding to the small fraction of 
quasiparticles that do relax to the bottom of the trap. 
These quasiparticles reside in the box until they are ther- 
mally excited by phonons. In the slow relaxation limit 
the activation exponent ([59|) can be reduced to 



1 u{5E) 



exp 



(65) 



[Rigorous evaluation produces a difference in the numer- 
ical factor here compared to Eq. (pT|) .] The tail of the 
distribution function describes the processes that 
are much slower than l/Fout, thus it must be retained 
despite its small amplitude, see also Eq. ([7T]) . 

The function F{t) defined in Eq. (|6Q|) can be evalu- 
ated using the small-z asymptote of rodd(^)7 see Eq. (|57|) . 
This approximation substantially simplifies F(t)^ 



F(t) 



dz 



tEo 



rr, 



■ exp 



out 



(66) 



Here we assumed that the main contribution to the F(t) 
comes from small-z region, z <C 5E/2T^ which limits the 

applicability of Eq. (j66|) to t <C T~^^ (|f )^^^- The asymp- 
totic expression for F{t) in Eq. (|64j) can be obtained using 
the saddle-point approximation 



F(t)'. 



.1/3 



exp 



2/3 



(67) 



The integral (|66|) can be also expressed in the analytic 
form in terms of the Meijer's G-function^^ . As one can 
see from Fig. [3 at low temperature T <^ SE there is a 
time window 



Tout ^ Tout \2Ty 



(68) 



in which the survival probability deviates from the expo- 
nentially decaying function. We assumed in Eq. (|68|) that 
the upper limit is more restrictive than t <C (Eout'?")^ 
so that r-dependent term in the exponent of Eq. (|67|) can 
be neglected. 

The fractional power 2/3 in Eq. (|67j) stems from the 
peculiarity of superconducting density of states at low 
energies. Assuming the quasiparticle distribution in the 
lead is given by Eq. ([3Q|) . every time a quasiparticle tun- 
nels into the box it may occupy a different energy level, 
which is reflected in the initial conditions, Eq. (|36|) . How- 
ever, due to the singularity of the escape rate T {Fk) 
at Fk ~ £^thd5 this results in a strong energy dependence 
of the dwell time of a quasiparticle. Therefore, averaging 
over many such events leads to the deviation of F{t) from 
the simple exponential function, as shown in Fig. H 



InF(t) 




FIG. 7: (color online). Deviation of F{t) (solid blue line) 
defined in Eq. (|6Q)) from the exponentially decaying function 
at Foutt > 1. (We assumed r = oo here.) 



At t > iw)^^'^ minimum of the exponent 
in (|66|) is beyond the limit of applicability of the small- 
z approximation for the rate Fout(^) given by Eq. ([57j) . 
and instead of Eq. ([66|) one should use Eq. ([6Q|) . Since at 
z ^ SF/2T the escape rate Fout(^) is a smooth function, 
F{t) decays exponentially. 



F{t)(x exp 



SE_ 
2T 



Tout (^min)^ 



(69) 



Here Fout(^min) - 

The lifetime distribution function A^odd(^) for the odd- 
charge state can be obtained from 5'odd(^) by substituting 
Eq. dm) into Eq. Under conditions of Eq. jM]) the 

lifetime distribution function Nodd{t) will deviate from 
the exponential distribution 



Nodd{t) 



24/3 



V3 (Foutt)' 



/3 



exp 



rout^ 



2/3N 



• (70) 



The average lifetime of the odd-charge state (Todd) 
the slow relaxation case is 



m 



(Todd) = 



-f 

Jo 



Sodd{t)dt 



1 



0rrFout7s 



7/ 



Despite the quasiparticle having small probability of re- 
laxing to the bottom of the trap, the main contribution 
to the average dwell time (Todd) is given by the tail of 
Sodd{t)' This is because once the quasiparticle is trapped 
in the CPB it spends there exponentially long time, see 
Eq. (|65|) . As expected (Todd) is the same for fast and 
slow relaxation cases since the average lifetime deter- 
mines the thermodynamic probability to occupy a given 
charge state. 



V. CHARGE NOISE. 

The complex statistics of capture and emission pro- 
cesses discussed in the preceding section also manifest 
itself in the spectral density of charge fluctuations of the 
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Cooper-pair box. In this section we study the charge 
noise power spectrum for "deep" (rPout <^ 1) and "shal- 
low" (rrout > 1) traps. 

The kinetic equations for occupational probabilities of 
odd- and even-charge state have the foim^ 

Peven {Ep,t)^^W (p, k) (Peven (^p , t) - Podd (^fe , t) ) = 0, 
k 

Podd {Ek,t)^^W (p, k) (Podd {Ek , t) - Peven {Ep ,t)) = 



(72) 



Here P^^di^k.t) = p^dd(^fe)^odd(t) with cTodd{t) = 
"^Zk Podd{Ek,t) , and the quasiparticle transition rate 
W{p,k) is defined in Eq. (|6]). Assuming that the lead 
is a heat bath of quasiparticles we can write even-charge 
occupational probability as Peven{Ep^t) = f{Ep)(Jeven{t) 
with f{Ep) being the distribution function of the quasi- 
particles in the lead, and creven(^) = "^p Peven{Ep^t) be- 
ing the occupational probability of the even state. This 
allows us to reduce Eqs. ([72]) to 



<^even 

k,p 

Podd{Ek,t) ^^W{p, k) {Podd{Ek,t)-f{Ep) it)) = 



(73) 



One can see that Eqs. ([75|) satisfy the normahzation con- 
dition: 



(74) 



The stationary occupational probabilities aeven and ^odd 
are given by the Gibbs equilibrium state. Assuming that 
f{Ep) is given by Eq. (|3Q|) . we obtain 



1 



^'"even — 



l + n^pHexp(^ 



CTodd = 1 — CTeven- (75) 



Here n^p is the quasiparticle density in the lead, see 
Eq. ([3ID, and is the volume of the CPB. 

The fiuctuations around this equilibrium state can be 
taken into account within the Boltzmann-Langevin ap- 
proach, which assumes that the occupational probabili- 
ties fiuctuate around the stationary solution ([75]) due to 
the randomness of the tunneling and scattering events as 
well as partial occupations of the quasiparticle states^^. 
The kinetic equations for the charge fiuctuations can 
be derived by properly varying Eqs. ([73]) and adding 
Langevin sources corresponding to the relevant random 



events^—, 
d 



+ Fin j (5aeven (t) = ^ ^ (p, ^) ^^odd {Ek , t) + ^ {t) , 

^ k,p p 

Wip, k)f{E^)5a,.Ut) HI (t) +Cr {t). (76) 



Here the relation (5creven(^) = — (^crodd(^) was taken into 
account. The Langevin sources ^J(/.)(^) and ^^^(t) cor- 
respond to quasiparticle tunneling from/to the state 
\p) I \k) through the junction, and inelastic electron- 
phonon scattering, respectively. [Note that ip (^) = 

— Xlfe ^fc (^) S/c S^(^) ~ ^-1 These random processes 
are considered to be Poissonian with the following corre- 
lation functions 

p 

= 2S{t — t')Sk,k'Tout{Ek) f {Ek — SE)ae^en, 



{ek'mi'it')) = s{t-t')^^^^^ ( Sk,k' ^^^^^^^ 



C^odd 

(77) 

The latter expression is consistent with the collision in- 
tegral in the relaxation time approximation and conser- 
vation of the probability crodd(^) by the electron-phonon 
scattering^^^^'^. 

The spectral density of charge fiuctuations in the CPB 
is defined as 

Sqiij) = 2e^(Jaeven(^)(^Creven(-^)), (78) 

and can be obtained from Eqs. ([76]) and ([77|) . The so- 
lution of the second equation of the system ([76]) in fre- 
quency domain is 

rp fij^ \ ^ont{Ek)f{Ek-SE) - ^pQ^^jEk) 

0Podd[Ek,U)) = : \ 1 OCTeveni^) 



-iuO^Tont{Ek)^\ 

g(^) + rM 



1 • 



-iuj + Tont{Ek) + - 
Substituting this expression into the equation for 

(5<Jeven(^) WC find 



^ -iw + rout(-Bfc) + - 



(79) 
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where the function C{(jo) is given by 



-iuo + - 

Y,f{Ek-SE) 



-iuj + Tont{Ek) + 7 



(81) 



Finahy, using Eqs. ([78|) and ([80|) we can find the corre- 
lation function ((5<Jeven(^)^^creven(— ^)), and obtain charge 
noise power spectrum Sq{(jo)^ 



2e2 



C{uj)C{-uj) 



E 



{lo^ {u^n, (-u;)) +ro,t(Jgfc)rout {Ew ) {it (^)g^ (-^)) 
(-iw + rout(^fe) + i) (iw + rout(^fc') + i) 



(82) 



Upon substituting correlation functions ([77|) into Eq. ([82|) 
the general solution for 5'q(cc;) can be obtained (after 
cumbersome but straightforward calculations, see the 
Appendix). Rather than going through the full deriva- 
tion, we study here Sq{(jo) in the limiting cases rPout ^ 1 
and rPout ^ 1, which can be derived from Eqs. ([77|) . 
(ED, and (|82l). 

We first consider fast relaxation limit rPout <C 1. In 
this case one can neglect the second term in the numer- 
ator of Eq. ([82]) . For cjr <C 1 one can simplify Eqs. ([8T]) 
and ([82]) further. After straightforward manipulations 
one finds that the leading contribution to the noise is 
given by Eq. ([23]) with the rate 



1 

Teff 



7f ■ 



(83) 



which includes all processes changing the population 
^even- The first term in Eq. ([83|) corresponds to the rate 



of thermal activation of the quasiparticles by phonons 
from the bottom of the trap to the lead, see Eq. (p!6|) : the 
second term is the rate of quasiparticle tunneling from 
the lead to the box given by Eq. (pT]) . also cf. Eq. ([24|) . 

In the opposite limit rFout ^ 1, the charge noise power 
spectrum Sq{(jo) can be roughly viewed as the superpo- 
sition of two Lorentzians, see Fig. [51 The first one corre- 
sponds to the processes involving quasiparticle thermal- 
ization and activation by phonons, and is dominant at 
low frequencies. The second (quasi) Lorentzian describes 
the fast processes (cj^Fout) associated with the escape 
of unequilibrated quasiparticles from the box. 

At low frequencies uo <C cc^cr, see Fig.[8l the noise power 
spectrum is well approximated by the Lorentzian func- 
tion. This can be obtained by neglecting the first term 
in the numerator of Eq. ([82]) . and keeping the leading 
terms in 1/rFout and oj/Fout in Eqs. ([8T]) and ([82]) . see 
the Appendix. After straightforward manipulations one 
finds 



Sq{uj) ^ 4e^Crodd (1 - CFo. 



1-D 



dd) 



Teff 



1 + C (Cjreff)2 + 1 

The constants C and D here are defined as 



• (84) 



C 



aA To, 



and D 



1 p{5E) 



exp 



^q(O) 
10-1 



10- 



10" 



10- 



10- 



10- 



10- 



1 



FIG. 8: (color online). Spectral density of charge fluctuations 
generated by quasiparticle capture and emission processes in 
the Cooper-pair box for the slow relaxation case (rFout = 
10^). Here <^cr ~ v^Fout/r is a crossover frequency between 
two different regimes governed by Eqs. (|84p and (|87p . 



The width of the Forentzian 

1^1 Tin 

Teff rFin + V^F, 



is given by 



out 



-+7s 



out 



out 



(86) 



The first term here corresponds to the transitions from 
even- to odd-charge state involving the relaxation of a 
quasiparticle to the bottom of the trap. [cf. Eq. ([25]) : 
difference in the numerical coefficients comes from rigor- 
ous solution of Eqs. ([77|), ([81|) and ([82|).] It is determined 
by the quasiparticle relaxation rate 1/r times the por- 
tion of the time the unequilibrated quasiparticle spends 
in the box. The second term in Eq. ([86]) describes the 
transitions odd to even state involving the escapes of a 
thermalized quasiparticle from the CPB by phonon acti- 
vation. It is proportional to the phonon-assisted quasi- 
particle escape rate from the box to the lead 7^ times 
the probability to ffnd an empty trap upon the escape of 
the thermalized quasiparticle. This probability is deter- 
mined by the portion of the time the trap spends in the 
even state upon the escape of the thermalized quasipar- 
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0.01 - \ 

0.1 0.2 0.5 1 2 5 10 Cj/Fout 

FIG. 9: (color online). The deviations of the charge noise 
power spectrum Sq{co) from the Lorentzian function at high 
frequencies u ^ Tout- Blue solid line corresponds to Sq{(jj) 
given by Eq. (|90)) . red dashed line is the normalized Lorentzian 
function with the width Tout- 

tide, and is determined by the fast processes involving 
Tout and Tin- 

At high frequencies, co ^ cjcr, the dominant is the first 
term in the numerator of Eq. (|82]) . Then, in the leading 
order in 1/ rFout the power spectrum becomes 

o (^Ar^^^^ CZi{u;)aeven . . 

(l + CZ2M)V(r^) (CZiH)' 

Here the sums over momentum k in Eq. ([82]) are replaced 
with the integrals (T ^ Sb). In terms of the dimension- 
less variable z ([52]) these integrals are denoted as Zi{uj) 
and Z2{u;) [see the Appendix], 




As shown in Fig. [51 in the frequency window Ucr <C Cc; <C 
Tout the noise power Sq{uj) becomes flat with the ampli- 
tude 

At higher frequencies u > Fout and C <C 1 the noise 
power spectrum ([87]) can be approximated by 

SQ{u;)^^CZi{uj)ae^en (90) 

I out 

with Ziiuo) given by Eq. ([88]) . At these frequencies the 
charge noise power spectrum Sq{uo) describes charge fluc- 
tuations due to the tunneling of the unequilibrated quasi- 
particles from the box to the lead. By taking a Fourier 
transform of Eq. ([90]) , one can notice that the noise power 



spectrum in time domain has the same functional form as 
Fit) of Eq. ([67]) . Therefore, charge noise power spectrum 
also reveals the deviations from the conventional Poisson 
statistics due to the singularity of the quasiparticle den- 
sity of states at low energies. The deviations of the charge 
noise power spectrum ([90]) from the Lorentzian function 
at cc;~Fout are illustrated in Fig. ([9|). At higher frequen- 
cies ix) ^ Fout charge noise power spectrum Sq{uj) decays 
as 1/cj^ see Eq. 1^. 



VI. CONCLUSIONS. 

In this work we studied the kinetics of the quasiparti- 
cle trapping and releasing in the mesoscopic supercon- 
ducting island (Cooper-pair box). We found the life- 
time distribution of even- and odd-charge states of the 
Cooper-pair box. We also calculate charge noise power 
spectrum generated by quasiparticle capture and emis- 
sion processes. 

The lifetime of the even-charge state is an exponen- 
tially distributed random variable corresponding to the 
homogenous Poisson process. However, the lifetime dis- 
tribution of the odd-charge state may deviate from the 
exponential one. The deviations come from two sources 
- the peculiarity of the quasiparticle density of states in 
a superconductor, and the possibility of quasiparticle en- 
ergy relaxation via phonon emission. The odd-charge- 
state lifetime distribution function depends on the ratio 
of the escape rate of the unequilibrated quasiparticle from 
the box Fout and quasiparticle energy relaxation rate 1/r. 

The conventional Poissonian statistics for both quasi- 
particle entrances to and exits from the Cooper-pair box 
would lead to a Lorentzian spectral density Sq{uj) of CPB 
charge fluctuations^^. The interplay of tunneling and re- 
laxation rates in the exit events may result in deviations 
from the Lorentzian function. In the case of slow quasi- 
particle thermalization rate compared to the quasiparti- 
cle tunneling out rate Fout, the function Sq{uj) roughly 
can be viewed as a superposition of two Lorentzians. The 
width of the first one is controlled by the processes in- 
volving quasiparticle thermalization and activation by 
phonons, while the width of the broader one is of the 
order of the escape rate Fout- 
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APPENDIX A: POWER SPECTRUM OF 
CHARGE NOISE. 

Combining Eqs. ([Ilj), ([80]), and ([Hj) we obtain the ex- 
pression for the charge noise power spectrum, 



Sq{u;) 



4e" 



(Al) 



{Lo^ + ^)rout{Ek)f{Ek-SE)ae.^n + rout{Ekrpldd(Ek)^ , a. 



odd 



(u;2 + (r„„t(£;fc) + l/r)2) 
Here the product w) is given by 



r,utiEk)ptddiEk) 



^ -iw + rout(i^fc) + 7 



C{l0)C{-L0) =L0^ 



p^dd(^^fe)rout(i^fc) , V- fiEk-SE)rout{Ek) 



-(rout(^fc) + l/^)' a;2 + (rout(^fe) + l/T)'^ 



' 1 P'odd (-gfc )rout {Ek ) (Fout (^^fc ) + 1/r) ^ f{Ek- 5E)ront {Ek){iV-' + l/TH Tout {Ek)/T) 



u;2 + (r,„t(i3fe) + l/r)^ 



(A2) 



Equation (|Aip can be simphfied in the thermodynamic Here C and D are given by Eqs. ([85]) . Substituting 
hmit by introducing functions Z\{bj) and Z^iyj) Eqs. (|A2p - ()A4p into Eq. (|Aip one can obtain the general 



^ ^ cj^ + (rout(£^/c) + l/r) 

and 

P^dd(^fc)rout(^fc)^ (^4^ 

^ V^' + (rout(i?fc) + l/r)' 

I 



expression for Sq{uj) 



4e2 


(rout) ^(rrout) 


'J- out 'J- out 




Tout / a; \^ 


i-^ZiH+cz2H 


2 


"^^2H+ 




2 



The functions Zi{lj) and ^2(61;) can be written in the The dimensionless variable z here is defined in Eq. ([52]) . 
form of the dimensionless integrals Assuming that at low temperature the main contribution 



to the integrals (|A6p and ()A7p comes from the small z 
region, z <C 5E/2T^ one 
K^io ''^c.2 + (roddW + l/r)^' ^""^^ using Eq. ([511) to obtain 



ry ( X Todd ("^ , e ^i/(2;)rodd(^) . region, z <C 5E/2T, one can simplify Zi(cj) and Z2(cj) 



and 



<'(ml« a;2 + (r.dd(j) + l/Tf 'Wo (a;/r.dd)!j + (l + v/J/7r„ddf 
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and Finally, by taking the appropriate limits in Eq. ()A5p 
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